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Q^l Abstract 

Mh. We calculate the leading twist contribution to near-forward proton-proton 

r~| ! (and proton-antiproton) elastic scattering with large momentum transfer, in 

the multiple scattering (Landshoff) mechanism. The amplitude in the near- 
forward region is dominated by singlet exchange for all three valence quark- 
quark scatterings. We assume the existence of a hard singlet quark-quark 
cd I amplitude, which we estimate to be 0{al/t). For a three-quark state whose 

transverse size is less than 1/Aqcd; Sudakov resummation accounts for both 
approximate da^^/dt ~ t~^ at moderate t, and da^^/dt ~ t~^^ at larger t. 
The transition from approximate t~^ to t~^^ behavior is strongly correlated 
with the transverse size of the valence three-quark state in the proton. 
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1 Introduction 

Perturbative QCD (pQCD) offers a formalism in which to study high energy ex- 
clusive hadronic processes fll, H, H- Factorization theorems make possible the sep- 
aration of long distance effects, describing the formation of color singlet hadronic 
bound states, from the short distance interactions of quarks and gluons. Once long 
distance interactions have been modeled by hadronic wave functions, the remaining 
hard process can be calculated perturbatively. At leading power, the relevant par- 
tons are the valence quarks that carry the momenta of the incoming and outgoing 
hadrons. 

Fixed-angle hadron-hadron elastic scattering has been discussed in this formal- 
ism, for both a single hard scattering [|l], |^, ^ and for multiple, independent ('Land- 
shoff'), scatterings of the valence quarks @, 0. To calculate the amplitude, however, 
it is necessary to incorporate systematically the effect of soft gluon exchange, as de- 
scribed in refs. [^ |^. This extended factorization can also be useful in the treatment 
of hadronic form factors p, |^ and proton-antiproton annihilation |Ty] . The result- 
ing suppression of spatially separated partonic scatterings accounts for approximate 
"quark counting" power behavior |I[] for fixed angle hadron-hadron elastic scatter- 
ing. In this paper, we will show that the improved perturbative formalism, including 
the transverse structure of the proton, can naturally account for the approximate 
t^® behavior of dcr^^/dt at small angles. In addition, we suggest a relatively sharp 
transition to dimensional counting behavior, dcr^^/dt ~ t~^°, as t increases. The po- 
sition of this transition is directly related to the transverse size of the quark valence 
state of the proton. 

Large angle scattering [9 = 0{'k/2)) involves two scales, the momentum transfer 
—t ~ s, and the scale of hadronic masses. In contrast, near-forward scattering 
{9 — > 0) is a three scale problem, involving the center of mass energy, y^, the 
momentum transfer —t and the scale of hadronic masses. Correspondingly, the 
near-forward amplitude at each order has both "Sudakov" double logarithms in 



momentum transfer and "Regge" logarithms of the energy |Tl|, |T2|. In a recent 



paper, we have studied the relationship between these two phenomena 



In this paper, we use the results of Refs. |T^ and to investigate the interplay 
between soft gluon corrections and the tranverse structure of the proton in a range of 
momentum transfers, 3 GeV^ ^ — t ~ 40 GeV^ . For the description of the hadronic 
bound state we use quark distribution amplitudes derived from QCD sum rules 



1^ , pig , |T6|, |l3|- To leading power in ag, the three valence quarks from each proton 
scatter independently. Of course, at the lower end of the t-range considered here not 
all individual q-q scatterings can be really hard. Nevertheless, the Sudakov-improved 
pQCD approach remains well defined, and affords a model for the amplitude over 
the entire range. 

In the s — > oo, t fixed, limit q-q elastic scattering is dominated by singlet ex- 
change in the t-channel. The suppression of octet exchange may be described as 



gluon reggeization |T2[. In [^ we have studied this process as the 9^0 limit of 
fixed angle q-q scattering. There, we rederived the reggeization of octet exchange 
from the energy- dependence of its Sudakov suppression, while the corresponding 
factor for the singlet is s-independent in the leading logarithmic approximation. Of 
course, even with singlet exchange, the elastic scattering of colored objects is IR di- 



vergent. In ITBI we have suggested ways of separating dimensionally regularized IR 



singular contributions from the hard part of singlet exchange. The renormalization 
group can be used to resum the leading and non-leading ln(t) dependence of the q-q 
amplitude, employing one-loop anomalous dimensions P, |l^, |l^ . 

In the following, we shall treat hard singlet exchange as an a^ process [|^ only 



(denoted H'^iQm) below), where m = 1,2,3 labels the hard q-q scattering of mo- 
mentum transfer Qm- Once the q-q scatterings are embedded in the full hadronic 
amplitude, IR divergences cancel, due to the color singlet nature of the external 
hadrons. The IR cutoff (denoted 1/6 below), although arbitrary when q-q scatter- 
ing is considered in isolation [jl3| , has a physical meaning in the case of proton-proton 



scattering, as the transverse separation between independent hard scatterings. 

The expiclit form of the hard singlet to lowest order, a^, is ambiguous, because 
the singlet amplitude is already IR divergent to this order. As pointed out in ref. 



T3[| , different IR subtractions procedures yield different expressions for H^. The IR 
subtractions that define the hard singlet to lowest order can be fixed only by their 
effect on the amplitude at higher orders in as- Without facing this difficult question 
here, we estimate the lowest order hard singlet in section 3 below. 

Because, as noted above, the transverse separation between hard scatterings 
plays the role of an IR cutoff, the transverse structure of the proton valence state 
wave function is an essential feature of our calculation. The perturbative (small b) 
behavior of the wave function is available from 0], but for the range of t considered 
here we will find it necessary to include as well a "nonperturbative" or "intrinsic" 
transverse structure as well as the familiar light cone wave functions. The impor- 



tance of transverse structure has been stressed in pi] , |22| , p3 |. 

Our discussion is organized as follows. In section 2 we describe in some detail the 
factorized form of the p-p amplitude in the multiple scattering scenario for fixed 
scattering angle 6, stressing the roles of hadronic wave functions and flavor flows. In 
section 3, we take this amplitude to the forward region, considering the contribution 
from singlet exchange in all three q-q scatterings. In section 4, we calculate the 
lowest order three-singlet amplitude. This calculation determines the combination 
of the quark distribution amplitudes that appears in the p-p amplitude. Leading and 
non- leading ln(t) corrections from all orders are then included through soft gluon 
resummation factors. In section 5, we define the kinematic cut-offs and the choice 
of scales that turn the closed form of the p-p amplitude into a numerically stable 
four dimensional integral, which we compute for a range of momentum transfers. 
We examine the rate of convergence of the amplitude in impact parameter space 
and its sensitivity to scale choices, and study the effect of the transverse structure 
of the valence three-quark state. We summarize our results in the final section. 



2 The proton-proton amplitude in the multiple 
scattering scenario 

We consider the leading twist contribution of the valence three-quark Fock state for 
each proton to the proton-proton elastic scattering amplitude. The p-p amplitude, 
A, is shown in fig. 1 for multiple scattering [^]. The valence quark momenta are 
parameterized in the CM frame as 



ki = ^J^Xiv'^ + i^ivt + K > (1) 

where i is the proton label, the Vi are lightlike vectors along the proton directions, 
with Pi = Js/2 Vi, and the v'^ are lightlike vectors in the directions opposite to Vi 
(fj ■ fj' = 1). kf = k^j_ denotes the components of ki in the plane perpendicular to 
the spacelike component of Vi. Using the notation Vij = Vi ■ Vj, we parameterize the 
invariants of the scattering as 

t = -SVi3 = -SV2A, U = -Sf 14 = -SV23 ■ (2) 

Since each proton contains two w-quarks and one (i-quark in the valence state, 
there is more than one flavor-distinguishable way of connecting the quarks to the 
three hard scatterings. We label the inequivalent flavor flows as follows. We call H^ 
the hard scattering blob in which the (i-quark from proton 2 participates, and H^ the 
hard scattering in which the d-quark from proton 1 participates, if H^ 7^ H^. Then 
we obtain the three inequivalent flavor flows shown in fig. 1, labeled by / = 1, 2, 3. 
Notice that for / = 1, H^ and H^ are flavor indistinguishable. 

For identical quarks, both t and u gluon exchange channels are available. For 
/ = 2 the gluonic exchange channels in the u-d hard scatterings [H"^ and H^) are 
fixed to be t, and for / = 3 they are fixed to be u. Table 1 summarizes the channel 
configurations of the hard scatterings for each flavor flow. 



2.1 The factorized amplitude 

According to ref. P], the leading twist elastic p-p amplitude may be written in the 
following factorized form, 



(27r 



|8 /•! 



3 

i i 



(3) 



where aj,aj are Dirac and SU(A^c) indices respectively, hi are proton helicities, 6 is 
the scattering angle in the CM frame, t] is the spacelike unit vector perpendicular 
to the scattering plane and i? is a remainder suppressed by a power of ^i. The 
y's in eq. (0) represent the valence structure of the external protons, the H"^^s 
represent the three hard scatterings and U summarizes the effects of soft gluons, in 
terms of two transverse momentum components normal to the scattering plane, i 
and i. From the kinematics of q-q scattering when all the valence quarks are on 
shell, Xj = ^1 > 0, Xj = ^2 > in eq. (|l[) for every i = 1, 2, 3, 4. The antisymmetric 
color structure of each wave function has been separated from the Y's and absorbed 
into the soft-gluon function U. U has remaining free color indices that are summed 



against those of the hard scatterings, H'^, as in Ref. ||I3[ . 

The hard scatterings depend only on the large components of the quark momenta. 
For this reason we may integrate the wave functions over the small components Ki 
and ki. The transverse momentum integrals in eq. (^, however, are linked by delta 
functions. The transverse delta functions in the factorized form of the amplitude 
reflect the physical picture of p-p scattering, as it procedes in the multiple scatter- 
ing mechanism. The incoming and the outgoing proton wave functions are highly 



Lorentz contracted along their respective light cone directions. Therefore, thought 
of as discs, they intesect along the direction perpendicular to the scattering plane 
0. The hard scatterings may take place anywhere along this line of intersection. 
Thus, in configuration space the wave functions are unlocalized in the normal to 
the scattering plane. When their separations are large, soft gluons moving in the 
normal direction can resolve the color of the quarks. These effects are summarized 
in U. The amplitude is conveniently represented by rewriting eq. (^) in "impact 
parameter" space in the transverse components. 



(27r 



6 ,1 



A{s,t,h,) = {-^)^ f d^,d^2 f dhdh 
sm ti Jo J 



3 



E E i;vL(^i'^2,&i,&2;P.,/i.) 



/ 1 ai,l^iy-fi 



;f^{a.b.cj(&l,&2,/i) n ^K.ajteo^i,/^), (4) 



m=l 



where /i is the renormalization scale and where Y, defined by 

t%,(6,6,&i,&2;P.,/^.) = /d/.,d«:,|d2k,d2k,e-(^-'^^+b-''^) 

^^aift7,(^«5 ^«5 -Pi) hi) , (5) 

is Y integrated over the k's and Fourier transformed with respect to its transverse 
momentum components, in terms of the vectors bj = bit], with rj the normal to the 
scattering plane. This is the factorized expression we shall employ below. 
Let us now discuss the various functions in eq. (|]) in more detail. 

2.2 Wave functions on and off the light cone 

Before assigning quarks to the independent scatterings we review the description of 
the (uud) Fock state. We omit the flavor flow index (/) in Y until we reinsert it in 
A, eq. (Bl). It is convenient to discuss the wave functions first in momentum space. 



As usual, the three-quark wave function is the Fourier transformation of the 
three-quark operator |^. For a proton moving in the -|- direction, with energy E, 
we have 



^ al3^\'^Ui ; "'M2 ) -^ ) '^5 A''/ 



NJ J (27r)4 (27r)4 



x(0|TK(yi)<(2/2)rf^(0)]|P,/i)6,,,. (6) 

In the hniit yf — > yfv^, the three operators approach the hght cone, and Y is related 



to a normal light-cone wave function ||2J]. /i is the factorization scale. As mentioned 
above, the F's are defined as color-invariant functions (appropriate ordered expo- 
nentials between the quark fields are suppressed). 

We have already seen that the fields in eq. @ are not on the light cone in general. 
Nevertheless, the ratio of transverse to longitudinal momenta of the quarks will, in 
all the wave functions discussed below, vanish in the high-energy limit. As a result, 
the leading twist contribution to the amplitude, eq. (^, may still be given by a 
spinor basis on the light cone. 

Permutation symmetry between the two M-quarks and the requirement of total 
isospin 1/2 imply that Ya^^ can be expanded directly in terms of products of spinors 
with definite helicity, 

M^X = (^1^2^3/2)-'/' Uc,{ku„+) up{ku„ -) d.,{kd, +) , 
A^g]^ = (E1E2E3/2)-'/' u^{ku„ -) up{ku„+) d^ikd, +) , 
M'^X = -(^1^2^3/2)-'/' u^{k,„+) up{k,„+) d^{kd, -) , (7) 

as 

Yap, = ^TJI^Mf^) [*123-MS, + ^213^S, + (^132 + Vl/231) -MSJ , (8) 

where the subscripts on the scalar function \1' refer to the order of momentum 
arguments, with ks = kd, for example, \l'(A;i, A;2, ^3) = ^123- In (|^), Ei, E2 and E3 
are the energies of the first and second w-quarks and the d-quark respectively. For 



a proton with (— ) helicity, it is only necessary to flip the spinor helicities in eq. (0). 
Here and below, spinors are normalized according to 

N{P, h)-f''N{P, h) = 2P^ . (9) 

A commonly-used alternate basis for Y is described in the Appendix. 

We are now ready to discuss the transverse structure of the three quark state, 
since this is crucial to the factorized amplitude as given in eq. (^). At high en- 
ergy, the dominant radiative corrections come from soft gluons, whose couplings are 
indpendent of the quarks' spin states p. The effects of these gluons on the wave 
functions factorize from the light cone quark distribution amplitude (referred to as 
LCDA below) and exponentiate in the impact parameter (6) space respresentation. 
That is, Y, eq. (|^), the Fourier transform of the the function \l'(fc„i, ku2, kd, /i) may 
be written as 



Yap-^{ii,i2,hi,h2]Pi,hi) — 2i/48A^ ^f^il^) 



^123-MS, 



+P213A^S^ + (Pl32 + P231) A^S^ 



(10) 



where 



^(^i,6,6,^i,&2;/i) = </'(^i,6,6;yu)^(6,6,&i,&2) 

= fdKif d/€2 / d^kie-^'^i-^^^i-^^) / d^kse-*''^-^^*^-''^) 

x^(A;i,A;2,A;3;/i), (11) 

with ?/^ is a dimensionless function of the transverse separations. To be specific, we 
denote by bm the positions of the hard scatterings H"^ along the 77 direction and by 
bm their mutual transverse separations defined as 

61 = 62 - ^3 , &2 = &i - ^3 , &3 = ^2 - &i • (12) 

Another connection between the function V and the LCDA is [^ 

^(^1, 6, ^3, h = b2 = 0; ^) = 0(^1, 6, 6; f^) + (a,(/i)) , (13) 



that is, neglecting perturbative corrections for large /i, 

^^(6,6,0,0) = !. (14) 

Alternately, in momentum space, \1/ is related to the LCDA by 

J dKi J dK2^{ku^,ku^,kd;n) =0(^i,6,6;At)^(^i,6,6,k„i,k„2), (15) 



with 



j d^kid^ks ^(^1, 6, 6, ki, k2; /x) = 1 . (16) 



At very high energy, the behavior of V in the 6j is a computable exponential 
0, which vanishes when any of the 6, equal 1/A. The perturbative transverse wave 
function, then, is of finite size. Rather than give its explicit form here, we have 
simply absorbed it into the resummed soft-gluon factor specified in eq. (PB| ) below. 
In the following, we shall assume that this has been done, and consider the function 
tp{C,i,bi) as a measure of the remaining "nonperturbative" transverse structure of 
the proton. As we shall see in Sect. 5 below, this structure has strong influence on 
near-forward elastic scattering. In a sense, we may think of ip as the source of the 
"intrinsic" transverse momentum of valence quarks in the proton. We should keep 
in mind, however, that the distinction between perturbative and nonperturbative, 
or pro cess- dependent and intrinsic, is somewhat arbitrary, and must depend, in 
particular, on the order to which the perturbation expansion has been carried out. 

As a last topic in this subsection, we briefly review the evolution of the light 
cone distributions. In the standard analysis 0, the LCDA is expanded in the basis 
of the Appel polynomials Aj as 



,2/A2\"/37 



/Inu^/A^ 
0(xi, X2, xs; jj,) = 0as(a;i, X2, xs) Y^ Nj ( ^— J ajAj{xi,X2, x^) . (17) 

0as (3^1,3^2,3^3) = 120x1X2X3 is the asymptotic limit of 0(/i) for /i — i> 00. /io ~ IGeV, 
(3i = (ll/3)A^c — (2/3)n/ is the one-loop coefficient of the QCD /5-function and 



A = Aqcd is the QCD scale parameter. The normahzation coefficients AT,-, defined 
by 

Nj j dxidx2da;35(l - Xi - X2 - Xs) (t)as{x)A^j{x) = 1 , (18) 

and the anomalous dimensions bj//3i are the same for all models. Model dependence 
of comes through the coefficients aj. For a list of the constants and functions 



Nj , bj , aj , Aj{x) see table 2 and [1^, ^. The dimensional parameter /Ar(/i) is given 
by 



/^(/i) = /7v(/io) ( ^z%l^ ' /^(/"o) = (5.2 ± 0.3) X lO-^GeV^ . (19) 



'ln/iVA2\ W 

Both and /at depend mildly on /i due to QCD evolution. 

The quark distribution amplitudes identified above contain much of the soft 
physics of the problem. However, information about soft gluon exchange is also 
contained in the perturbative resummation factors. As mentioned above, these will 
include the perturbative b dependence of the wave functions p. To organize this 
dependence, we must analyze color fiow in hard and soft exchange. 

2.3 Color flow, hard and soft exchange 

According to the conventions of fig. 1, ki,ki,Pi — ki — ki are the momenta entering {i = 
1,2) or leaving (z = 3,4) the blobs H^,H^ and H^ respectively. The identification 
of these momenta with {ku^)i, {ku2)ii {kd)ii is obtained once the fiavor fiow / is fixed. 
Similarly, the third index 7 of the proton wave function, eq. (§), is identified with the 
Dirac index (3i or 7, in eq. (|^) according to the hard scattering if™, m = 2, 3, in which 
the d-quark of the corresponding proton participates. Each hard-scattering function 
if" consists of lines with virtuality ^^t, with ^^m '"^ = 1,2, 3, the longitudinal 
momentum fraction carried by the quarks in scattering m (they are all the same 
for a given m). 

The color fiow basis suitable for our treatment is the octet-singlet basis, 

10 



(CadjjjQ.} — (,-'mja3ai(,-'mja4a2 ) lCsj|^.j — 00103^0204 • (.-^Uj 

For near-forward scattering the anomalous dimension matrix associated with H"^ 
becomes diagonal in this basis []TB|, ITBI, M . More detail will be given in the following 



section. 

The color tensor U includes soft radiative corrections that cannot be absorbed 
into the proton wave functions. As such, its perturbative expansion begins at zeroth 
order in the couphng, 

4 

t/{a,b,c0(&l,&2,/i) = l[^a,b,c, + ^K) " (21) 

i=l 

At lowest order, U consists simply of the e color tensors associated with the external 
protons. In ref. the one-loop radiative corrections for fixed angle p-p scattering 
were analyzed, and the corresponding renormalization group equations identified. 
The IR safety of the p-p process was shown to arise from the cancellation of IR 
divergences among graphs of the form shown in fig. 2. The underlying physical 
principle is that there is an effective IR cut-off 1/R for the transverse momentum 
of the exchanged gluons, with R the distance between adjacent if^'s. Below this 
cut-off, long wavelength gluons cannot couple to the color singlet hadronic state. 

Leading logarithms in s//i^, t//i^ and bffi"^ may be resummed into exponential 
facotrs, exp(— s/), one for each hard scattering function ifj". As above, / = s, adj 
labels the color flow of H"^ in the basis (|20|) . Beyond one loop in the corresponding 
anomalous dimensions, the exponentials for each hard scattering no longer factorize 
from each other in general. Once leading logarithms have been organized in this 
fashion, the renormalization scale /i in H"^ is replaced by the hard-scattering scale 



Qm — ^mC2V —t , (22) 

where C2 is a parameter chosen to optimize higher order corrections. Similarly, in 
U , /i may be replaced by a scale of the order of the l/6j's. In the following, we shall 
neglect higher orders in as(l/6j), and retain only the first term in the expansion of 
f/, eq. dgg). 

11 



The exponents that are generated by the one-loop radiative corrections to the 
hard scattering H"^, with color flow /, are given by [^ 



^ J \Simi "m) 



Cl 



C3 



g„ln 



Xn 



Hm ~r C,n 



+ C2 



Hm 1 I Hm 

m 



Cn 



An 



2(ln(2g„) + 1) - 2^(ln(2C™) + 1) - \n\2qm) + ln'(2C„ 



c 



A \ Cm / 



with 



qm = ln(Q„/A) = \n{UC2V^/A) , 
The positive constants Ci, C2, C3, /5i and /52 are 



ln(|6^|A). 



(23) 



(24) 



Cl : 

C3 ^ 

/52 



8CV 
34 

y 



11 2 

A = y AT, - -ny 

-20 



^^^^ (18 ~ ^^^V ~ 9^'^''^ ^ ^f/?i(7e - In 2) 



(25) 



N'. - {^^K + ^C,) T,n, , 



For A'"c = 3 and normalization Tp = 1/2, we have Cp = 4/3 and we take Uf = 3. 
P2 is the two-loop coefficient of the QCD /5-function and •je is the Euler constant. 
Of particular interest is the last term B^^ in eq. (0), which contains the color flow 
information of the hard scattering, 



5, = A, + 2Cf In 



327i3-r 



- N,Cf 



(26) 



Here, A^ is the eigenvalue of the anomalous dimension matrix of H^ along the color 
flow direction /, and the term —NcCp comes from the anomalous dimensions of the 
participating quarks, 47^ = —NcCfoIs/t^ in the axial gauge. The Xj will be specified 
in the next section. 
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To summarize our results so far, the p-p elastic amplitude may be expressed in 
terms of integrals in 6-space in the general form 

(27r)6 ^1 



A{s,t,h,) = (-^)L^ / d^id^ /dfeidfea 
sm 6 Jo J 



3 



xE E >^„Vft,,(ei,6,&i,&2;P.,^. 






y<U{a.^b,^^}{bi, as{l/bi)) {cjj^^ (c/J^^ (c/3)^^ 
xexp(-S7^(Qi,6i) -s/2 ((51,62) - si.^{Qi,h 

X n(^/™)K}te'^-«^(Q™))- (27) 

m=l 

3 The quark-quark hard amphtude in the near- 
forward region 

In the kinematic region of near forward scattering with large momentum transfer, 

s > -t > A^ , (28) 

the leading power contribution to the p-p elastic amplitude comes from the "ttf'gluon 
exchange channel for / = 1 and / = 2 (see table 1). The information about color 
flow in each hard scattering if" enters through the eigenvalue A^ in the B^ term 



of the exponent, eq. (0). In |T^ we have calculated the one-loop octet and singlet 
eigenvalues. The result is 

Ai = Aadj = N, In ^—j - — + As , As = 2Cp, (29) 

giving 

5adj -B, = N, In (^^) + 0(ln° s) . (30) 

Consequently, the contribution to the p-p amplitude from octet exchange in a single 
hard scattering is suppressed relative to singlet exchange by a factor 



p^' 
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with Qm, Cm defined in eq. (p^. We note that the A/ in eq. ( p9D are gauge dependent, 
but that the full result in eq. ( p3D is gauge invariant 0. Due to eq. (0), we expect 
that the dominant contribution to A in the near-forward region comes from singlet 
exchange in all three hard scatterings, which we denote by ^3s. 

For the octet, if™[j begins at single gluon exchange, but the lowest order con- 
tribution to the singlet, H^ is at the one-loop level, and the contributing on-shell 



graphs are IR divergent. The IR finite hard part of H^ can be put in the form |13 

r'-MQlUlt) = ^^^'(^-) In (^) , (32) 

where 

N^ -1 

is the color coefficient coming from the projection of the lowest order color structures 
along the t-channel color singlet, and where we define 

p parameterizes the ambiguity in the IR subtraction to lowest order. We expect an 
optimal p to be determined by an examination of higher orders. 

Inspection of eqs. (p2D , for Qm, and (^), for r^, shows that an "MS" choice for 
the parameter C2, 

C? = -, ^ (35) 

gives a hard singlet that is proportional to p. We take for C2 this MS choice for the 
rest of this paper and denote 

r™(«.(C),dt) = -. ^^^;ff-V . (36) 

Sm'' 

With the choice (^) for C2, the Bg coefficient of the Sudakov exponent, eq. (p6|), 
using eq. (|^) for As, becomes 

B, = 2Cf{ie + In tt) - N,Cf ■ (37) 
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4 Calculation of the p-p amplitude 

At this stage, all the necessary ingredients are available for the calculation of ^33. 
First we calculate the lowest order (a^) contribution to the hard-scattering ampli- 
tude. The calculation is straightforward using eqs. (0), (||), and (|). To present the 
result in a compact form we define 

Tl23 = 2 [-^132+^231], (38) 



for each external wave function, with V defined in eq. (|ri|) . The result in the impact 

parameter representation, eq. (p7D, is found from 

2 _ 

Yl Yl ^ifft7.(^i' ^2, hi, 62; Pi, hi)JJ{aMc,}(bi, a.(l A)) (c/J„^ (c/J^^ {cj,)^^ 






m=l 



-2 ]2 
'I23J 



= -r WC<5,,,35,,,, /^(/i)^{[n'23 + ^213 + 4ri? 

+ [n'23 + ^213 + 4Ti223] [n'32 + ^Il2 + 4T1I2] } . (39) 

Notice that the first term inside the braces, from / = 1, is symmetric under 1 ^^ 2, 
because for / = 1, H^ and H^ represent the same flow of flavor. 

The color factor C is determined by the normalization N^l of the F's, eq. (^, the 
eikonal tensor U, eq. (0), and the color structure of the iJ^'s, eqs. (|3^) , (|33|) . It is 
(73/4 \ f N"^ -l\^ 1 

^ " (i\gj4 (^^^'^'''^'7 *^^''*{'*»} '^^"''i^'} *^^''*{^'} " (^ 4iv2 j ^i\g)2 ' '^^'^^ 

yielding C = 2/3^ for A^^^ = 3. 

Including the soft-gluon factors for singlet exchange and using 
sin^ e = -4t/s{l + 0{t/s)) we obtain 

AUs,t) = -jrT-2l deid^/ d6ld62[T3s](^)7^(el,6,&l,&2;/i) 
(47rj^ t Jo J-oo 

X exp[-Ss(Qi, |6i|) - s,(Q2, I&2I) - Ss(Q3, l&i " ^2!)] , (41) 
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with the lowest order contribution [Tss]^^-* given by 

= ~^(27r) -2-27" T -^^-^ pdh^h-A^hi- (42) 






The second hne in the above equation follows from eq. (0), which defines the hard 
scatterings, and we denote by Wm the hard scale at which each hard scattering H^ 
is evaluated. The wave function combination TZ is 

+ \Pl2M + 'PllM + ^'^UMVlz2{l^) + ^Il2(/^) + 4r,232(/x)]} . (43) 

In summary, the three singlet amplitude is given by eqs. (^TD-(^3D, with soft gluon 
resummation exponents for singlet exchange given in eqs. (ESI) and (P7|). These 



results are the same for proton-antiproton near-forward hard elastic scattering, since 
the corresponding amplitude is obtainable from s ^-^ u crossing of p-p scattering, 
and in the forward direction u^ k. s^. 

5 Numerical results 



The experimental study PB| , ^ of proton-proton near-forward elastic scattering 
reveals two characteristic features in the momentum transfer range of interest here. 
The first is the independence of the differential cross section dcr^^/dt from the energy, 
for an energy range i/i ^ 25 GeV up to 62 GeV. The second is the power behavior 
da^^/dt ~ t~^ . The lowest order multiple scattering model [^] predicts precisely 
this scaling and power behavior. Landshoff and Donnachie [^ were able to fit the 
experimental data in the range 3 GeV^ ^ — t ~ 15 GeV^ by 

^ = Ct-^ , C = 0.09 mb(GeV)^ . (44) 

at 

First, we examine the factorized amplitude eq. ( pTj) in this region without transverse 

structure for the wave functions beyond the soft gluon exponentials exp[— Sg], that 
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is, with ilj{^i,bi) = 1 in eq. (pJ^. Then we shall include a model for nonperturbative 
transverse structure and study its influence. 

5.1 Perturbative transverse structure 

Upon retaining only the LCDA's in the proton wave- function (see eq. ([T3|)), the 
function TZ eq. (^31) becomes 



^ = /^(/^)^{[0L3(/^)+0213(/^)+4Ti223(/^)]' 

+ [0L3(/") + 02i3(/^) + '^T^23if^)M,,{fi) + <Pl,M + 4T^,,{fi)]} , (45) 

with Ti23 = (l/2)(032l + 0132)- 

Energy- independence of the cross section follows from eqs. (|^), (^2[) and the 
relation 

dt IGtts^ 
in which the 1/s^ cancels with the s^ in |^p from eq. (|41|). Also, the singlet eigen- 



value As, appearing in the exponents Ss, eqs. (^), (poD, is energy independent, unhke 
the octet, Aadj- Exact s-independence is not anticipated to persist at very high en- 
ergies, where the energy dependence of the singlet exchange hard scatterings should 
become important |T^. 

The t~^-dependence is characteristic of the multiple scattering model if we as- 
sume that the integral over the suppression factors exp(— Ss) yields negligible t- 
dependence. Indeed, then Ass ~ s/t"^ and, from eq. (^BD, da^^/dt ~ t^^. This power 
dependence cannot persist at high enough t, where the strong Sudakov suppression 
in the exponents leads to the onset of dimensional counting behavior of the form 
^3s ~ s/t^ and hence to daPP/dt ~ t-^°. 

In the numerical study that follows we plot the quantity A^~^|^3s| t^/s as a 
function of —t. The dimensionless normalization coefficient A^ is 

A^ = (27r)^(47r)^--gC = 6.895 x 10^ . (47) 

Pi 
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It has been defined without the factor p , and we use the one-loop running couphng 
asiw"^) = 47r/(/5i ln(ti;^/A^)) in [Tss]'-^-*, eq. (^). The numerical results for the 
amplitude Ass are in units of \p^\. The straight line labeled 'DL' in the figures 



corresponds to the fit of eq. (0). In order to compute ^3s from eq. ( ^ll) we need to 
fix the 61, 62 integration limits and the scales fi,Wi,W2,W3 . This we do as follows. 
i) Lower hmits for 61, 62- The exponents Sg have a root at Q^ = l/|&m| (see eq. 
(^)). We refer to the region < \bm\ < 1/Qm as the small-6 region. Formally in 
this region the factor exp[— Sg] in eq. (^Tj) enhances instead of suppressing. We take 
its value in the small-6 region to be 1, i.e., 

SsiQm, \bm\) ^ S,(Q^, K\)e{Q^ - l/\bm\) (48) 

and integrate down to \bm\ ^0 0. 

ii) Upper limits for hi, 62- \bm\ are strictly limited to be less than 1/A, as past 
this value the Sudakov exponents diverge. We impose upper cut-off be in the 61, 62 
integrals and examine the behavior of Ass as b^ -^ 1/A. 

in) For the scales fi , Wm we make the following choices. The factorization scale 
/i, separating wave functions from the remainder of the amplitude, (eqs. ([T7|), (|19|)), 
is taken equal to the minimum average transverse momentum of the valence quarks, 
corresponding to the maximum transverse separation between the hard scatterings, 

u = min \ -^, -^, -^ ^ } . (49) 

\\bi\'h\'\b^-b,\j 

In the soft quark region, ^m — > 0, the hard scales Wm shift from Q^ to p. as 

Wm = raax{Q^,fi} . (50) 



In addition to the conditions (^), (|5(]|), we impose a lower bound /imin for all the 
scales 

fI,Wi,W2,W3> /imin, (51) 

and examine the sensitivity of the amplitude as /imin — *■ A. 
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The result depends, in principle, on the values of /imin and A, the model for 
(f){xi] fi) and the value of the undetermined parameter p in eq. (0). In this numerical 
study we take \p\ = 1, but it should be noted that the results for the amplitude scale 
with IpI^. We start with /imin = 1 GeV and the LCDA model of Ref. [l^ and look at 



the rate of convergence of Ass as b^ -^ 1/A, as well as its sensitivity to the choice of 
the scale parameter A. The results are shown in figs. 3, 4. Note that for A = 0.1 GeV, 
Aimin = lOA corresponds to 6 = 1/yUmin = 0.1/A and for A = 0.2 GeV, yUmin = 5A 
corresponds to 6 = 1/yUmin = 0.2/A. As the figures show, up to these values for be 
the Sudakov suppression has hardly started affecting the amplitude. Therefore, in 
the region of the impact parameter space where Sudakov effects become important 
the coupling constant is effectively fixed at as(/^min) ^^^ ^'^ condition (|5lD. 

We tentatively divide the t-range into three regions corresponding to small (a), 
intermediate (b), and large (c), absolute slopes. We observe approximate linear 
dependence of ln(|.43s| f^/s) on In (t) in each of the above regions. Table 3 lists the 
slopes of the integral for selected values of be- The slopes on the logarithmic plot 
correspond to 't'^'op'" deviation of the amplitude from the t~^ behavior. The onset of 
dimensional counting in region (c) is apparent, where strong Sudakov suppression 
stem the growth of the amplitude. Note also the relatively small deviations from 
the t~^ behavior in region (a), less pronounced in region (b), where the transition 
to dimensional counting occurs. 

The dependence of the amplitude on A is also of interest. We observe that, for 
this choice of /imin and LCDA model, A^s decreases with increasing A for b^ ~ 0.3/A 
while it increases past this value of be- This is due to the competition between two 
counteracting effects. The first is the decrease of the impact parameter space volume 
with increasing A due to the Jacobian factor 1/A^ in 

be _ _ 1 rbcA _ _ 

dbidb2 = -T^ / db[db'^ , (52) 

A^ JO 

noting that the integrand is a function of bcA only. The second is the increase of 

the coupling constant with A. From A = 0.1 GeV to 0.2 GeV, as(/^min) increases by 

a factor ~ 1.4. This gives a substantial effect because the lowest order contribution 
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to the amplitude is 0{a^). In the case of the proton form factor 0, where the 
lowest order contribution is (9(a^), the shrinking of the impact parameter space is 
always stronger than the increase of as with increasing A, even when /imin -^ A. 
Proton-proton near- forward scattering is much more sensitive to changes in a^, even 
though the Sudakov suppression here is stronger relative to the form factor, due to 
the presence of the positive eigenvalue As = 2Cf (eq. (p9D ) in the exponent Sg. 

Decreasing the lower bound /Xmin down to 0.7 GeV, we recompute the amplitude 
for the same combination of parameters as before. The results are shown in figs. 5, 
6. Again here, for be ~ 0.3/A the coupling is essentially fixed at asifJ'mm)- The main 
feature is the marked increase of the amplitude (notice the scale shift on the ^3s 
axis). For A = 0.1 GeV, A^s increases by a factor ~3.0 in region (a) and by a factor 
~2.3 in region (c). For A = 0.2 GeV the increase is by a factor ~5.5 in (a) and 
~3.4 in (c). This shows that the calculation is quite sensitive to the lower bound 
/imin- Although Sudakov suppression always makes the integral finite, removing the 
condition ( [51| ) ( equivalent to letting /imin -^ A ) leads to an increase by orders of 
magnitude. 

For fixed t, the relative logarithmic distance between the be curves remains almost 
constant with varying /Xmin at fixed A. This is a safe diagnostic that the dependence 
of the amplitude on /imm is due to the coupling constant, since for the largest part 
of the impact parameter space «« = as{n'^in) o^ lii~^(/imin/A). The drift of Ass to 
higher values with decreasing /^min does not affect the slope of the be curves in region 
(c). Dimensional counting is bound to set in for high enough t. But the transition 
region from approximate t~^ to dimensional counting behavior shifts towards smaller 
t with decreasing /imin- This is not due to any change in Sudakov effects, but to the 
fact that the smaller t region receives larger corrections from ^^(/imin)- 

The general characteristics of the ^3s behavior described above are the same 
for all LCDA models, and the following conclusion can be drawn. Sudakov sup- 
pression is not strong enough by itself to stabilize the pQCD calculation of the p-p 
near-forward amplitude against contributions from soft regions. This is not to be 
interpreted as a failure, but rather as an indication of the importance of the trans- 
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verse structure of the three-quark Fock state niodehng the proton. It is apparent 
that in the t-region where t~^ behavior is observed, the Sudakov suppression sets 
in at unphysical scales, larger than the transverse size of the proton. t~^ behavior 
for the amplitude is due to "geometrical" scattering, to which the entire size of the 
proton can contribute, and the transition to dimensional counting occurs only once 
the Sudakov suppression is strong enough to confine the independent scatterings to 
regions smaller than the average transverse size of the three quark state. In this 
way we can assign physical meaning to the scale /Xmin, as being of the order of the 
average transverse momentum of the valence quarks inside the proton. Following 
this reasoning, we proceed by including in the calculation the transverse structure 
of the proton. 



5.2 Intrinsic transverse structure 

Returning to the factorized form of the hadronic wave function, eq. 

we consider a transverse wave function of the form P, ^ ^, |^, [l^ 

ij{xi,ki) = Nexp 



3 1^2 



i=l -^^ 



(53) 



(54) 



where N = Ri/{7T'^xiX2Xs) and R^g is a mean impact parameter. The impact 
parameter representation of ip{xi, kj) is defined in eq. ([TT|) . The result is 



ip{xi,bi 



exp 



4^1, 



(X1X263 + xiXsbl + X2X361) 



(55) 



The value of A^ is fixed by the requirement ^(6j = 0) = 1 (eq. (|T^)). In this paper 
we do not advocate any particular choice of '?/'(xj,kj), but use (S3) because it has 
received much previous attention. 
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The assignment of quarks to hard scatterings according to the flavor flow / does 
not affect ilj{xi,bi), because it is totally symmetric in its indices. The expression of 
the amplitude, eq. (|41|), can be readily obtained in this case from TZ, eq. (|5), which 
now reads 

+ [0?23(/^) + 0213 (/^) + 4Ti^23(/^)][0?32(/") + 0312 (/") + ^Tl^2{^^)]} . (56) 

The parameter R^g is a measure of the transverse size of the three-quark Fock 



state. Brodsky, Huang, Lepage and Mackenzie [£T| have presented arguments on how 
to estimate it. Here, we keep it as a free phenomenological parameter and examine 
the behavior of the p-p amplitude within a range of values of R^g ~ 1/fcrms- A crucial 
point is that Rsg < 1/A. Since ip{xi,hi) contains the transverse size information we 
set the maximum of the h integrals to h^ = 1/A. 

The only, yet conceptually important, modification in the choice of scales /i, Wm 
is in condition ([51|), where /imin is replaced by l/Rsg, 



1 , , 

/i , wi , w;2 , W3 > — — , (57) 

since l/R^g is a typical transverse momentum. 

The numerical results for ^3s are shown in figs. 7-10 for the various LCDA models 
and for fixed A = 0.20 GeV. The numbers on the R^g curves are the slopes in their 
respective linear regions. 

Before drawing conclusions from the numerical results, we list some of the ap- 
proximations made in the course of this calculation, (i) we have integrated only 
along the direction perpendicular to the scattering plane, where the two Lorentz 
contracted disks intersect. This picture becomes less clear-cut as 6* — > and con- 
tributions from the integration over the azimuthal angle in the impact parameter 
space will have to be considered; (m) we have assumed that the gaussian, eq. (p^D, 
correctly describes the transverse momentum distribution; {in) we have neglected 
all hadronic mass scales (proton and constituent quark masses) that would normally 

22 



show up in the transverse wave function ip; (iv) we have neglected the transverse 
momenta in the hard q-q scatterings Hm [§|, by assuming that their contribution 
in the soft quark region can be captured by taking as{w^) = as{l/Rlq). We sug- 
gest, however, that these and other shortcomings, such as the neglect of higher-order 
corrections in U and the if™, do not affect the main features observed in figs. 7-10 
and shared by all LCDA models. 

In figs. 7-10, we see: (z) approximate, but not exact, t~^, behavior for the 
amplitude for moderately large t, where we find ^3s ~ s/t^'^. We note the fit 
daPP/dt ~ t-^-^ of the FNAL data, suggested by the authors of ref. ^. (ii) A 
transition from Landshoff to dimensional counting behavior, whose position is a 
sensitive function of the transverse size R^q. Since the experiment is decisive about 
the existence of Landshoff behavior at least up to —t = 15 GeV^, we can infer that 
our calculations imply a small transverse size for the valence state, i?3g ~ 0.3 fm 



2T|. A much larger value of Rsq would exhibit t ^° behavior in the cross section at 



lower t. 

Our results so far are encouraging, but we recognize that the energies and mo- 
mentum transfers involved are not so high that the perturbative picture is applicable 
without reservation. In particular, we would like to get an idea of the importance 
of soft quarks, and consequently low momentum transfers to our results. When one 
of the scales C,i becomes small, the corresponding momentum transfer decreases as 
C,f, and at some point the transverse size of the "hard" scattering may be as large 
as the separation between hard scatterings. 

Thus, as an additional test of the independent scattering method, we have recal- 
culated the amplitude with a cutoff .^j > A on soft interactions. Fig. 11 shows the 
effect of demanding that all three of the C,i be larger than A. The A = 0.0 curve is 
the reference, and corresponds to the calculations in the preceding figures. The con- 
tributions of soft quarks are substantial, but not overwhelming. Demanding that 
the softest quark have l/20th of each proton's momentum leaves more than half 
the amplitude, while demanding 1/lOth reduces the amplitude by a factor of about 
three. We should keep in mind, however, that at collider energies, even a quark with 
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^ = 0.05 is quite energetic. Of course, at such scales, the quark-quark momentum 
transfer is no longer very large, and a treatment that includes the transverse mo- 
menta of the wave functions as part of the hard scattering ^ is probably necessary. 
Fig. 12 shows the same calculation, but now with the restriction imposed only on 
the two more energetic quarks. Here we note that most of the time the middle 
quark has at least one-tenth of the energy, and it has more than l/5th about one- 
third of the time. This shows, we believe, that the independent scattering picture, in 
which there are more than one, and sometimes three, independent and well-localized 
scatterings, describes the basic picture of proton-proton elastic scattering near the 
forward direction. 

6 Conclusions 

We have described proton-proton near-forward elastic scattering with an improved 
pQCD factorization formalism, which describes both short-distance t~^° and Land- 
shoff t~® power behavior for moderate momentum transfers with s ^ — t. The 
presence of t~^ behavior in the cross section implies the relevance of a new scale, 
of the order of the average transverse momentum of the constituent quarks. The 
inclusion of a nonperturbative, or intrinsic, transverse size somewhat smaller than 
1/A seems to be required for the description valence state. Indeed, were the proton 
three-quark much larger than 1/A, Sudakov effects would produce t~^° behavior in 
the kinematic region where t~^ has been observed. All of these results require color 
singlet exchange for the independent quark-quark hard scatterings. 

We note that it is most natural for us to compare the calculated amplitude to 
experiment, rather than the cross section directly, since in elastic scattering, it is the 
amphtude and not the cross section that is the object of computation. Of course, 
the range of uncertainty in the cross section is much larger than in the amplitude. 
Nevertheless, it is quite striking that for |p| = 1 in eq. (|3^ the normalization of the 
amplitude works so well, considering the high powers involved; see eqs. (|19]), (^Tj). 

The transition from LandshofF to short-distance power behavior is strongly cor- 
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related with the transverse size of the three-quark state modehng the proton. We 
beheve that further understanding and technical improvements on the theoretical 
level combined with measurements of the position and the width of the transition re- 
gion could enable us to determine the average transverse size and estimate the form 
of the transverse structure from the experimental data. Evidently, near-forward 
elastic scattering involves spatially-separated hard scatterings of colored objects. It 
is clearly of interest to investigate a role for such phenomena in diffractive jet and 
other semi-inclusive cross sections. 

Despite the numerical significance of soft-quark contributions, the model de- 
scribed above displays a satisfying self-consistency. In fact, we know of no other 
model for p - p elastic scattering that describes the t~^ behavior of the elastic scat- 
tering cross section. We also think it more than ever of interest to determine the 
normalization and s dependence of the singlet scattering amplitude. 
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A Matrix wave functions 

As an alternative to eq. (0) , Y^is^ may be expanded in terms of spin matrix tensors 
E^'^K ci= 1,2,3 pj 



Y^p, = J2 Erp,{P, h)Y^''\K,,K,, k,- /i) , (58) 

for incoming protons (i = 1, 2) and in terms of a conjugate basis 

Y^p, = Y.Efp,iP. h)Y^''\K„K„ k,; /i) , (59) 



for outgoing protons [i = 3,4). The E^s and -E"s are given by 



Ei5^ = E-''\i^,C)^pN, , E%^ = -E-'l\C''^,iUN, , (60) 

£7 = ^E-'^\^,uV^CU{Yi,N), , EiX 



Eg, = zE-'/\a,,v''C)^p{j^j,N), , E^X = -^E-'/\C-'a^,v')^p{N^,r) 



We have defined these spin structures to be dimensionless. Here, C is the charge 
conjugation matrix, while 

(^i^v = ^[1,.,'yu], 75 = ^e^^KA7''7''7V , (61) 

with eoi23 = 1 and N{P, h) the massless proton spinor, normalized, as above, to 

N{P, h)-i^'N{P, h) = 2P^ . (62) 

The real functions F'-'^-' in this basis are conveniently normalized in momentum 
space as 

Y^''\ku,,k^,,kd■,^l) = 21/44JV | /^(^)^^''Hfe»i. fe»2, kd, /i) , (63) 



where, in a variation of the notation of ref. |]14[ for instance, 
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/Ar(/i) is a standard normalization constant, specified below. The primes indicate 
that these wave functions still depend on all four components fcf and have non-zero 
mass dimension. The various symmetries of the wave function require that V, A' 
and T' may be themselves obtained from a the single function \E'(fcu^, ku2, kd', fi), as 

V'{kui,ku2,kd;fi) = l/2[^ {kui, ku2, kd] fi) + '^{ku2,kui,kd] ft)], 
A [ku-^, ku2i kd] /i) = ^ 1 2y^ \ku2i km, kd', jJ.) — M/^fc^jj, ku2, kd', l^)\ , 
T'{ku„ku2,kd',iJ^) = l/2[^{km,kd,ku2',lJ^) + '^iku2,kd,km',fi)]. (65) 

Eq. (0) enables us to relate V, A' and T' to dimensionless functions of the light-cone 



variables ^j only (the functions V, A and T of Ref. p^). 

Linear combinations of the spin structures E^'''' parameterize the on-shell states 
with definite quark helicity defined by eq. (|^). The correspondence for our normal- 
ization choice, and for a proton with (+) helicity, is 

^i^j, = E«, - Ei5^ = {E,E2Es/2)-'/' u^{km,+) up{ku2, -) d,{kd, +) , 
M%, = E^^l^ + Eg^ = {E,E,E,/2r'/^ u^{km, -) upikm,+) d,{kd, +) , 
M'-X = E% = -(^1^2^3/2)"'/' Ua{ku„+) up{ku2, +) d^{kd, -) , (66) 



where, as in Sec. 2, Ei ,E2 ,E3 are the energies of the first and second u-quarks and 
the d-quark respectively. 

With this form of the wave function, the calculation of the p-p amplitude involves 
Dirac traces over the spin structures E^'^\ which are rather lengthy but straighfor- 
ward. The following relations are helpful [r] = (+, — , — , — )): 

(N3rN,)(N,'y^N2) = -2s6h,hA.h, (67) 

Y^tr{[7^ , 7^]7"[7p , 7<7]7^} = ^/.[p'^?'^^' + V^iJp^S'il'' - VpipVa^uV"'^ ■ (68) 
The result, of course, is the same as in eqs. (|4T|), (|42|) and (|^) above. 
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Table 1. 
The gluonic exchange channels in the hard scatterings for each flavor flow. 



/ channels in H^H^H^ no. of channels 

1 ttt, ttw, tut, tuu, + t ^^ u 8 

2 ttt, utt 2 

3 tuu, uuu 2 



Table 2. 
The decomposition coefficients in terms of Appel polynomials, and the anomalous di- 
mensions of 0(xj; /i), for the models of Chernyak and Zhitnitsky (CZ) |[l^, Chernyak, 
Ogloblin and Zhitnitsky (COZ) [|^, King and Sachrajda (KS) [|16| and Gary and 
Stefanis (GS) P 



J 


h 


aj (CZ) 


aj (COZ) 


aj (KS) 


aj (GS) 


Ns 




A,{x) 








1.00 


1.00 


1.00 


1.00 


1 




1 


1 


20/9 


0.410 


0.350 


0.310 


0.391 


21/2 




Xi -X3 


2 


24/9 


-0.550 


-0.424 


-0.370 


-0.588 


7/2 




3X2-1 


3 


32/9 


0.357 


0.460 


0.630 


-0.749 


63/10 


3(xi- 


X3)2 + X2(5X2-3) 


4 


40/9 


-0.0122 


-0.00259 


0.00333 


0.0176 


567/2 


(l/3)( 


Xi -X3)(4X2- 1) 


5 


42/9 


0.00106 


0.0633 


0.0632 


0.574 


81/5 


7x2 - 5 + (14/3) (x? + xl + 3siX3) 



31 



Table 3. 
The slopes of the 6c-lines in figs. 3, 4 considered approximately straight in the 
t-regions (a), (b), (c). 

slopes for A = 0.1 GeV slopes for A = 0.2 GeV 

KA (a) (b) (c) (a) (b) (c) 



1.0 


-0.38 


-0.61 


-0.84 


-0.23 


-0.52 


-0.89 


0.8 


-0.28 


-0.52 


-0.81 


-0.17 


-0.42 


-0.79 


0.6 


-0.19 


-0.39 


-0.70 


-0.10 


-0.33 


-0.69 


0.4 


-0.09 


-0.27 


-0.56 


-0.02 


-0.23 


-0.58 
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Figure Captions 



1. Proton-proton elastic scattering in the multiple scattering scenario, and flavor 
routings for each flavor flow /. The dashed lines represent the (i-quarks. All 
momenta flow from left to right. 

2. Graphs whose IR divergences mutually cancel. 

3. The three-singlet amplitude |^3s|, modulo p^Ns/f^, as a function of —t for the 

= 1.0 



CZ model |jTj], without transverse structure and for A = 0.1 GeV, /in 
GeV. 



4. |^3s|, modulo p^Ns/t^, as a function of —t for the CZ model |T^, without 
transverse structure and for A = 0.2 GeV, yUmin = 1.0 GeV. 



5. \A3s\, modulo p'^Ns/t^, as a function of — t for the CZ model |T^, without 
transverse structure and for A = 0.1 GeV, /imin = 0.7 GeV. 



6. IvAssI, modulo p^Ns/t"^, as a function of —t for the CZ model |n], without 
transverse structure and for A = 0.2 GeV, yUmin = 0.7 GeV. 



7. |^3s|, modulo p^Ns/f^, as a function of —t for the CZ model IQ, with trans- 
verse structure and for A = 0.2 GeV. The ISR data are from ref. [^ at 



^/s = 52.8 GeV, and the FNAL data from ref. ^ at v^ = 27.4 GeV. The 
error bars include both statistical and normalization errors. 



8. |^3s|, modulo p^Ns/t*, as a function of —t for the COZ model [|T^], with 
transverse structure and for A = 0.2 GeV. 



9. |w4.3s|, modulo p^Ns/f^, as a function of —t for the KS model |jT6|, with trans- 
verse structure and for A = 0.2 GeV. 



10. |v43s|, modulo p^Ns/t^, as a function of —t for the GS model fl^, with trans- 
verse structure and for A = 0.2 GeV. 
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11. |.4.3s|, modulo p'^Ns/t'^, with lower cutoff A on the scale ,^mcd of each of the 
hard scatterings. 

12. |^3s|, modulo p^Ns/t"^, with lower cutoff A on the scale .^mcd of the intermediate 
hard scattering. 
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